Quantum Theory of particles and fields as an extension of a probabilistic variational 
approach to classical mechanics and classical field theory. I 
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A theoretical scheme, based on a probabilistic generalization of the Hamilton's principle, is elab- 
orated to obtain an unified description of more general dynamical behaviors determined both from 
a lagrangian function and by mechanisms not contemplated by this function. Within this scheme, 
quantum mechanics, classical field theory and a quantum theory for scalar fields are discussed. As a 
by-product of the probabilistic scheme for classical field theory, the equations of the De Donder-Weyl 
theory for multi-dimensional variational problems are recovered. 
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I. INTRODUCTION 

The main purpose of this paper is to explore the pos- 
sibility of reconsidering the problem of the relativistic 
quantum theory of fields from a different point of view. 
We will focus on standards fields with no reference to 
strings. 

Let us start from the following general remarks. Sup- 
pose to have a classical isolated physical system described 
by a set of real space-time functions (fields) q^{x){i = 
1,2, ..,n), X — {x'^ , , x'^ , x^) E M, the space-time man- 
ifold, with the metric tensor g^^, = diag{\,—l,—l,—l). 
If it is associated to the system a Lagrangian (density) 



C^C{q\x),d^q\x)) id. 



dxt^ 



,/. = 0,l,2,3) (1) 



then, as is well known, the field equations determining 
the functions q' (x) can be obtained through a variational 
principle, that is by requiring that the action 



Ar = 



Iv 



C{q'{x),df,q\x))dx [dx = dx°dx^dx'^dx^) 

, (2) 

be stationary under variations of q''{x) vanishing at 
the boundary of V, a domain of M. Then one obtains the 
partial differential (Euler-Lagrange) equations 



dC 



d{df_,q^{x)) dq'ix 



= 



(z = l,2,..,n) (3) 



In the case of a classical mechanical system with a 
finite number of degrees of freedom (m.s.f.), a useful 
and relevant tool for the study of its dynamical prop- 
erties is given by the Hamiltonian approach. There is 
essentially one Hamiltonian formulation associated to a 



'Electronic address: ' matteo.villaniQba.infn.Ttl 



given Lagrangian. This formulation has a fundamen- 
tal role in the transition from classical to quantum the- 
ory, if this transition is performed on the basis of the 
standard canonical quantization rules (s.c.q.r). But it is 
also very relevant and effective in conjunction with the 
Hamilton- Jacobi theory and the general aspects of the 
one-dimensional variational problem involved in the case 
of a m.s.f. 0,0,11. 

On the other hand, in the case of continuous systems 
described by fields, we do not have a univocal Hamil- 
tonian formulation associated to a given classical La- 
grangian. The multi-dimensional variational problem ^ 
is tackled essentially according two different points of 
view: one developed mainly in physics, the other in math- 
ematics. 

In physics, a continuous system is interpreted as a me- 
chanical system with an infinite number of degrees of free- 
dom (m.s.i) [1, Q. The spatial coordinates {x^,x'^,x^) 
have the role of labels of these infinite degrees of free- 
dom. The Hamiltonian formalism is developed on the 
basis of this mechanical transcription. Such a formula- 
tion has a crucial role in the transition to the quantum 
theory of a continuous system, since it allows to apply 
the s.c.q.r. of quantum mechanics. However it has a 
draw-back in the fundamental distinction between time 
and space variables. As is well known, this problem is 
bypassed in quantum theory, by working in the Heisem- 
berg or Interaction picture, or in the framework of the 
Feynman's functional approach. But, due to the original 
point of view, further problems appear, like the occur- 
rence of infinite quantities. In some cases these problems 
are bypassed through additional prescriptions, like the 
renormalization. In these cases, the resulting (effective) 
quantum field theory has been very successful. It has 
reached the stage of a paradigm, although accompanied 
by doubts about its completeness. It runs into difficulties 
in the case of the gravitational fields. 

At a classical level, a draw-back of the mechanical 
transcription of a continuous system is that it leads to 
an Hamilton- Jacobi equation formally of functional type, 
which frustrates its use. On the other hand, the math- 
ematical analysis of multi-dimensional variational prob- 
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lems does not make an essential use of the interpretation 
of a classical continuous system as a m.s.i. [Ull, In 
this analysis, the space-coordinates are treated on the 
same footing of the time: we have four parameters which 
together take the place that the time alone has in a classi- 
cal m.s.f. As a matter of fact we have several approaches 
to multidimensional variational problems P, , which 
share this point of view. However one relevant approach 
has been formulated, called the De Bonder- Weyl Theory 
[i1)[1,S01j which introduces and makes use of an "Hamil- 
tion" function. In this theory a covariant Hamiltonian 
formalism is elaborated Q and a generalized Hamilton- 
Jacobi equation is obtained, which is not a functional, 
but a partial differential equation. 

At the classical level, the De Dondcr-Weyl formalism 
appears more appropriate than the Hamiltonian formal- 
ism based on the mechanical transcription. However, if 
we consider the transition to a quantum theory, it is in 
conflict with the s.c.q.r. . In fact these rules require only 
one canonical momentum (density) associated to each 
field that is d£/doq^{x) (as happens in the the- 

ory based on the mechanical analogy), while in the De 
Donder-Weyl theory we have four momenta 



dC 



(a* = 0,1, 2,3; z= 1,2,.., n) (4) 



Several attempts have been made in order to overcome 
the above difficulty and to obtain a covariant canonical 
quantization [Tol [Til , and reference therein] . Unfortu- 
nately, in these attempts, there is some arbitrariness and 
some insubstantiality in the probabilistic interpretation. 

It seems that the mechanical point of view is necessary, 
if one insists on the s.c.q.r. as the unique tool for a 
quantum theory of fields. In the pas t there are been 
objections to this point of view As a matter of 

fact. Born [l^ tried to avoid the s.c.q.r. and to introduce 
an alternative quantization procedure based on a result of 
the multi-dimensional calculus of variations (the Hilbert 
independence theorem) p^ . But his attempt was not 
successful. 

In this paper we make again an attempt like that of 
Born, but firstly we elaborate an appropriate alternative 
theoretical scheme for the analysis of the transition from 
classical theory to quantum theory. Developing an ap- 
proach proposed in a previous work , we adopt at the 
start, that is already at the classical level, a probabilistic 
point of view. In this setting we make a systematic use 
of the variational approach, by introducing a generaliza- 
tion of the Hamilton's principle. We obtain in this way a 
scheme more fiexible and more directly open to possible 
extensions. 

The preliminary aspects of this scheme are discussed in 
Sectini in the context of classical mechanics. The tran- 
sition to quantum mechanics is analyzed in Sect lllll this 
transition appears more as a natural extension than a 
conceptual jump. A particular role is covered by time- 
reversal invariance and by the local or global character of 
the dynamical equations. We will deduce the request for 



a mechanical system to be natural and discuss a mecha- 
nism for the emergence of complex amplitudes. The case 
of a discrete random variable (spin) is also analyzed. The 
probabilistic point of view is applied to classical fields in 
Sect llVI We avoid the mechanical transcription, since it 
leads to an ill-defined probabilistic scheme. However we 
show that there exists a well defined probabilistic scheme, 
in which the usual deterministic classical field theory can 
be embedded. It is very interesting that, in this scheme, 
there appear, as a by-product, the equations of the De 
Donder-Weyl theory. We see that the two Hamiltonian 
formulations, which are "equivalent" from the determin- 
istic point of view (in the sense that they give the same 
equations have different implications if we adopt a 
probabilistic point of view. 

The transition to a quantum theory of fields is dis- 
cussed in Sect|Vl This transition, in which the s.c.q.r. 
are avoided, is an extension of the results of the sections 
miandllVl 



II. PROBABILISTIC VARIATIONAL 
APPROACH TO CLASSICAL MECHANICS 

We consider a conservative mechanical system with n 
degrees of freedom and configuration space R" , described 
by a Lagrangian L, i. e. a given smooth function of 2n 
independent variables 

L = L{q,w) [q^ (g\ g^ .., g"), w = (i«\ w^, .., w")) 

(5) 

The coordinates q give the configurations of the system. 
Let us fix a time t > and make in ([5]) the substitutions 



q ^ m , w ^ i{t) 
Then we can consider the action functional 

/= r Lm.mdt 

Jo 



(6) 



(7) 



According to the Hamilton's principle, the temporal evo- 
lution of our system is describes by its trajectories q(t), 
which are functions of t such that 



— L{q{t)+er^{t),q{t)+em)dt 



= 

for every smooth function rj{t) having support in [0, r]. 
As is well known we obtain in this way the Lagrange 
equations for q{t) 



(i = l,2,..,n) 

(9) 

Now, we adopt at the start a probabilistic point of view. 
At each instant t e [0,r], we consider the configuration 
q a random variable described by a probability density 
p{q,t). We limit ourselves only to this assumption, that 
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is we do not require a priori that the temporal evolu- 
tion is necessarily described by a stochastic process. We 
will take into account of the obvious fact that the total 
probability 



p{q, t)dq [dq ^ dq^.. dq^) 



does not depend of the time t, assuming that this con- 
servation law can be expressed in a local form 



dtp{q,t)+dkf{q,t)^0, 



(10) 



through a probability current density j{q,t) — {j^{q,t), 

We will consider j{q, t) as a further dynamical variable, 
which, together with p{q,t), characterizes the "state" of 
the system. As a matter of fact, ([TU|) takes the place or 
generalizes the link between q{t) and q{t) of the stan- 
dard theory based on the notion of trajectory. We plan 
to determine (adding initial conditions) p{q, t) and j{q, t) 
directly from the given Lagrangian ([S]), through an ap- 
propriate generalization of the Hamilton's principle [15| |. 
To this end, suppose that p{q,t) and j{q,t) describe an 
actual temporal evolution of our system and let us con- 
sider the family of subsets E — U x [to,ii]j where U is 
a compact domain of E" and 0<<o<^i^''"j such 
that p{q,t) > for (7,t G S. In each S let us make the 
substitution 



w 



(11) 



in the given Lagrangian ([5]). Due to the meaning of j(g, t) 
and p{q,t), this substitution is a generalization of 
Then let us consider the action functional 



dtdq 



L{q,^^)p{q,t) 



X{q,t){dtp{q,t)+dkf{q,t)) 



(12) 



where X{q,t) is a Lagrange multiplier field defined in E. 

We can now formulate a generalized Hamilton's prin- 
ciple, which replaces ([5]), as follows: 

"The actual functions p{q,t) and j(q,t) are such that the 
equation I JO)) holds in R" x [0, r], and in each E there 
exists a function X(q,t) such that 



d_ 



j + £i3o,P + £2Po,X 



(i = l,2) (13) 



for every smooth function jQ{q,t) for which jo{q,t) = 
for q G 9U, and every smooth function po{q,t) such that 
Poiq,to) = Pa{q,ti) = and pQ{q,t) ^ for q <E dV." 
By considering we deduce from this principle 



dkX{q,t) 



dL . i{q,t) 

Qyjk V 



p{q.t)) 



), (fc = l,2, ..,n), (14) 



while from §^ = we obtain 

f{q,t) dL , j{q,t) 



dtXiq,t) 



p{q,t) dw^ ' p{q,t) 



(15) 



in each E. We note that, in the derivation of p4p and 
the condition p{q, t) — for q € 9U has no role. 
However this condition has some implications which will 
be discussed in the following. 

We can solve ([T4l) with respect to j/p and substitute 
the result in (fTU)) and To this end, we assume that 

our Lagrangian satisfies the condition 



det{ 



iq,w))^0 



(16) 



We see that the variation of p in the generalized action 
principle leads directly to consider the classical Legendre 
transformation. If we introduce the Hamiltonian H 

H{q,p) =w''pk - L{q,w) , {p = {pi,p2, ..,Pn)) (17) 

where w — w{q,p) is the solution of the system of equa- 
tions 

dL 

Pt^ -Q—{q,w), (i = l,2, ..,n) (18) 

we obtain, from and P7|) . 
f dH 



p dpk 
Then ITSl) becomes 



{q,dX) (5A = 5iA,..,5„A) . (19) 



dtX + H{q,dX) = (<z,tGE) 



(20) 



In the following we will call global any property, equa- 
tion, function and so on, which is valid or defined on 
the whole M" x [0,t], more precisely for every G 
R" X [O, t]. Otherwise we will make use of the term local. 

In the evaluation of the global dynamical variables 
p{q, t) and j(g, i), to which we are interested, we are faced 
with the field variables A(g, t), which, in principle, can be 
local. This is related to the presence or absence of a set 
of points where p{q,t) — O. But, a priori, we have no 
knowledge of this set (which we will call N). Further- 
more the fields X{q,t) have to satisfy, in any case, some 
consistency conditions, since they are related, through 
(|19p . to the variables p and j. However in the equation 
pop we don't have any term depending on p (and then 
on the condition p{q,t) > O). Therefore, as a result of 
this decoupling of X{q,t) from p, the lack of knowledge 
of N has no infiuence. The consistency conditions can be 
solved automatically by promoting the X{q, t) to a global 
function S{q,t), thus satisfying the global equation 

dtS{q, t) + H{q, dSiq, t)) =0 (g, < G M" X [0, r]) (21) 

From we deduce that j/p can be promoted to a 
global variable determined by S{q, t). Then we can write 
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for p the global equation 

n TT 

dtp{q,t)+du{p{q,t) — {q,dS{q,t)))^Q (q, t £ M" x [0, r 

(22) 

The equation ([21]) is the well known Hamilton- Jacobi 
equation 

0,S,S11- Its decoupling from ^ IS respon- 
sible for the basic aspects of classical mechanics. In fact, 
the equation ([22]) describes a deterministic process. It ad- 
mits, as particulare solutions, probability densities p{q^ t) 
having the structure 

p{q,t)=5{q-q{t)) , (23) 

where q{t) satisfies the system of first order differential 
equations 

BIT 

q\t) = ^{q{t),dS{qlt),t)) (fc-1,2, ..,n) (24) 

For this g(t), one can deduce, from (pij) and the 
canonical Hamilton's equations and the Lagrange equa- 
tions We see that the type of randomness determined 
by a pure Lagrangian approach is essentially related to 
initial conditions in the configuration space. 

From the start we have avoided deliberately references 
to trajectories or else to ensemble of trajectories, with 
the aim to have a scheme more suitable for possible ex- 
tensions. If we make use of trajectories in the proba- 
bilistic approach, as happens if we assume initially that 
j{q,t) — p{q,t)v{q,t)) {v G R"),with v{q,t) indepen- 
dent of p and regular in R" x [0, r], (a relation which is 
deduced in the previous scheme), then the equation (PT|) 
and (|22p can be also obtained through alternative well 
known classical procedures fl^, [l3, ■ In this respect 
we note also that (|14l) and (|15p . with j/p substituted 
by the so called "geodesic fields" (or control fields) have 
been already obtained by Caratheodory in his approach 
to existence problems for the extremals of the classical 
action 

©fllilll. They have been called by him the 
fundamental equations of the Calculus of Variations J^] . 
These ingredients of the Caratheodory approach appear 
in a natural and unified way in the previous probabilistic 
scheme. 

But the essential advantage of the generalized action 
principle is that it provides a springboard for a suitable 
more general scheme. The aim is to have an unified de- 
scription of more general dynamical behaviors involving 
both deterministic processes given by a Lagrangian func- 
tion L and processes not contemplated by this function. 
As a matter of fact, here we will not be concerned with 
external (random) actions on our system, but with more 
general internal mechanisms which govern its dynamics. 

The possibility that there can be active additional 
mechanisms is strictly related to the structure of our gen- 
eralized action principle. In general we can implement 
the property that the total probability is independent 
of time, through a balance equation involving dtpiq, t) 
and several quantities connected with the dynamical pro- 
cesses taking place in our system. It can be expressed in 



a form like that of a master equation. So we can natu- 
rally extend the previous classical scheme by considering 
balance equations more general than the equation pU)) . 
Hence, our approach is based in general on two ingredi- 
ents, a Lagrangian function L and a master-like equation, 
which can involve additional processes that cannot be as- 
cribed to this function. 

However, it is useful to take into account of a further 
aspect of the generalized action principle, which allows 
to see classical mechanics from a new perspective. It 
is concerned with the non uniqueness of the Lagrangian 
function L. 

In the standard variational approach, based on the tra- 
jectories, the arbitrariness of L is expressed by the addi- 
tion of a total time derivative. We have a similar situa- 
tion in the case of the generalized action principle, which 
is strictly related to the validity of the conservation law 
(fTU)l . In fact, due to pH]) . the Lagrangian L and the 
Lagrangian L' given by 

L' = Liq, I) + i^j^dkxiq, t) + dtxiq, t) , (25) 

where x(9:^) is an arbitrary smooth function, arc equiv- 
alent. They give, through the action integral (fT2|) . the 
same equations of motion. Therefore the choice between 
L and L' is irrelevant. However we have another class 
of equivalent lagrangians, involving only our internal dy- 
namical variables p{q,t) and j{q,t). It can be easily ver- 
ified that, due to pI7|) . also the Lagrangian L given by 

L = L{ql) + d{p)^-^^^dkp{q.t) , (26) 

where d{p) is a given function of p, is equivalent to L. 
This result is strictly related to the condition that the 
variation of p(q, t) must vanish on the boundary of U. 
This condition appears superfluous at the classical level, 
if we take into account of the structure of the action func- 
tional p^ . However the equivalence of L and L allows 
us to interpret the scheme of classical mechanics from a 
new point of view. In fact, if we consider L as our La- 
grangian function, then we can say that at the classical 
level there are operating also virtual internal dynamical 
mechanisms, generically of diffusive type. Their manifes- 
tation is forbidden by the particular relationship between 
p and j, given by the balance equation (|10p . 

The fact that (fTU]) makes irrelevant the coupling be- 
tween i and dp appearing in (|26p . is related to a funda- 
mental invariance property in the dynamical behaviour of 
a mechanical system. Let us assume that the Lagrangian 
L satisfies the condition 

L[g,w)^L{g,-w) , (27) 

as happens in general for the kinetic energy term. Then, 
as is well known, we have that the law governing the 
evolution of our mechanical system, given by the classical 
scheme, is invariant under time-reversal. Formally, the 
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system of equations (|T0)) and (|20|) (with j given by ^19\\ ) 
is invariant under the (time-reversal) transformations 

t' = -t 



t 



q" = q' 



p{q,t) ^ p'{q,t')= p{q,t) 
X{q,t) ^ X'{q,t') = -X{q,t) 
q,teJ: q,t' eT.' = V X [-ti, -to] (28) 

Now the term j^dkP in (|26p ) is not invariant under the 
transformations (p8| V Its ineffectiveness in the classical 
scheme can then be ascribed to the action of the time- 
reversal invariance. 

III. QUANTUM MECHANICS 

We consider now a particular and relevant extension 
of the classical scheme. 

First of all, in the framework of the generalized action 
principle, we take L as the Lagrangian function of our 
system, that is we take into account of the virtual in- 
ternal processes which are present at the classical level. 
However we don't assume that dkj^ gives all the contri- 
bution to the time derivative of p{q, t). We can interpret 
the field j{q,t), which through j/p determines the value 
of L at the point q, as a quantity related in a broad 
sense to the transport or convective aspects of the dy- 
namics. Together with p, j{q,t) characterizes, as before, 
the "state" of our mechanical system. But now, in the 
balance equation, we take into account of an additional 
non convective contribution to dtp{q,t) by considering, 
besides j, a current density jdiq, t) having a given explicit 
dependence on p. Here we limit ourselves to a balance 
equation which can be still expressed in a local form. So 
we generalize the equation to 

dtp{q, t) + dkjHq, t) + dkjiiq, 0=0, (29) 

with id a certain local function of p 

jd{q,t) ^i{q,p{q,t),dip{q,t),didjnp{q,t), ...) . (30) 

As a matter of fact, through the current jd in the bal- 
ance equation we consider a real additional mechanism 
in the dynamical evolution, which is, in a broad sense, 
of diffusive type. However it makes real also the pro- 
cesses contemplated by L, which are virtual when jd = 0. 
Starting from L and the extended balance equation (f^ . 
the action integral in the generalized action principle be- 
comes 



dtdql \L(q, -) 
s 1 P 



d{p)^-dup\p{q,t) 
P 



+ \{dtp + duj'' + dki\q.p.dip,..))\ (31) 



It. [j, P, A] 



From (|3ip we deduce the generalization of (fTil) and , 
that is 

BT i 

dkX-dip)dkP^j^{q/-) (fc = l,2, ..,n), (32) 



+ -^iq,-)-L + dip)dkj' + {dkX)^ 
p ow'^ p op 



dk{{diX) 



d{dkp)' 



dkdi{{dmX) 



d{dkdip)' 



= 0, 
(33) 



in every E. 

By introducing the Hamiltonian H , we obtain the sys- 
tem of the two coupled equations 

dtX + H{q,dX~d{p)dp) + 

+ d{p)dk{p^^{q,dX-d{p)dp)) + 

+ (*A)--a(,aA)_) + 



-I- dkdl{{drnX) 



didkdip) 



)-... = (34) 



OH 

dtp + dk[ p-^iq,dX- d{p)dp)) 



(35) 



in every S. 

In general, in the system of equations ((34l) and (|35p . 
we lose time reversal invariance. However we impose to 
our extension of the classical scheme the explicit condi- 
tion that the invariance under the transformation P5|) 
is preserved, notwithstanding the presence of the current 
i{q, p, dip, ..) and of d{p)dp. As a matter of fact we have 
two types of diffusive processes, which can have a compet- 
itive or antagonistic role. As we will see, the requirement 
of time-reversal invariance selects an important class of 
mechanical systems. Moreover, it allows to deduce com- 
pletely, from the knowledge of L, the specific structure 
of jd- 

A condition for the invariance under the transforma- 
tions (pS)) is that d{p) must be necessarily different from 
zero when jd 7^ 0, as can be easily seen from the equations 
([55)1 . Moreover, by considering we deduce from 
that this invariance requires 



2 
dp 



^ ^ <^(p)dp) - ^ d{p)dp^ 



H{q, dX - d{p)dp) - H{q, -dX - d(p)dp) 
+ d{p)d, 
+ {dkX) 
+ dkdi{{d^X) 



dk[{diX) 



ddkpj 
- ... = 



d{dkdip) 



From the structure of the equation ([36|) and the re- 
quirement that it must be satisfied identically by A and p, 
we deduce that H{q, dX — d{p)dp) — H{q, ~dX — d{p)dp) 



(36) 
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must be linear in dX, or, in other words, that H{q,p) 
must be quadratic in p. So, the request that the system 
(p4| . ([35|l be invariant under the time-reversal transfor- 
mations (pS)) . forces us to limit ourselves to natural [2^ 
Lagrangian systems, that is to mechanical systems de- 
scribed by a Lagrangian of the type 

L^L{q,w)^ {q)w'^w^ - V{q) , (37) 

where mij{q){^ mji{q)) is a positive definite mass ma- 
trix. 

From we deduce also 



so that 



d{dkdip) d{dkdidjp) 



jd = dp). 







(38) 



Hereafter we will fix our attention on the Lagrangian 
([57|. We have 

H{q,p) = ]^m'^{q)p,p^+V{q) , (39) 

where m*-' is the inverse matrix of mij. The condition 
of time-reversal invariance becomes 



+ dkdjX 



pd{p)m'^^ - 



di 







(40) 



In order that (PD|) be satisfied identically by A and p, 
we must have a relationship between d(p) and the current 
«(g, /9, 9p). It is given by 

l\q,p,dp) ^m^={q)pd{p)d,p + l\q) (fc = l,2, ..,n) 

(41) 

It can be easily verified that (|^T|) makes also the equa- 
tion l|35p time-reversal invariant, provided dki^{q) = 0. 
So «*'((/) is irrelevant and can be completely neglected. 
We conclude that i'^ must be a diffusion current of the 
type 



-D''^q,p)d,p, 



(42) 



with a matrix diffusion coefficient determined by the 
mass matrix and by the " coupling" function d{p) (so that 
in general it depends on p), i. e. 



D'^^{q,p)^~pd{p)m'^^{q) 



(43) 



So, we have worked out a formulation of a class 
of dynamical behaviors, which generalizes the classical 
scheme, but preserves in any case the property of time- 
reversal invariance. This class is characterized by the 



action functional ([3T|) . the Lagrangian (|37|) and the re- 
lationship (|¥l]) (with i'^ — 0). The law governing these 
dynamical behaviors follows from and ([55]) . We have 



dtX + ^m'^ d.Xdj X + V+]^m^^ {pd^ [p)) 'd.pdjp 
- d,{m'^pd^p)d,p) ^0 , 



dtp + dXpm'WjX) = 



(44) 



(45) 



in every S where p > (( )' = ^())- 

It is useful to note that the system . has an 
Hamiltonian structure, with p(q, t) and X{q, t) playing the 
role of conjugate fields variables. From (|?T|) and ([5^ it 
follows easily that (|44l) and (|45| can be deduced directly 
from the standard variational principle applied to the 
action functional 



At. [A, p] = / dtdq 



pdtX - He 



(46) 



with the effective Hamiltonian density Tie given by 



Tie{p,dp,X,dX) = +p 
1 



im'='9fcAa,A + V 



+ ^pd^{p)m'''dkpdip. (47) 

(The variation of A must vanish on the boundary of S). 
A particular case of (ITfl) has been previously postulated 
in the Fenyes theory [2l| . 

Now, our problem is to determine the global field vari- 
ables p{q,t) and j{q,t), having at our disposal the two 
coupled local equations (|44|) . (|45|) and the local relation- 
ship ([5^ . From this we deduce, as before, that the local 
(in principle) fields X{q,t) must satisfy some consistency 
conditions. In this respect we observe that if pd^{p) is 
a smooth function of p regular at p = 0, there appears 
no effective restriction on the set of points in R" x [0, r] 
in which the equations (|44p and (|45|) can be considered 
valid. In this case, the Hamiltonian Tip is a smooth func- 
tion in every domain of M" x [0,t]. As a consequence, 
in this case, the A(q, t) can be promoted to a global field 
variable S{q,t), as happens in the classical scheme. But 
now the equation for S{q,t) involves also p{q,t). Such a 
situation will be called the extended classical domain. We 
can make a further step and go beyond the extended clas- 
sical domain, by considering the situation in which pd^{p) 
is singular at p = 0. In this case the above straightfor- 
ward globalization does not apply. However we can ask if 
this singularity can be absorbed in a canonical transfor- 
mation of the field variables p{q,t) and X(q,t), allowing, 
in terms of the new variables, a straight-forward global 
description of dynamical behaviors not contemplated by 
the extended classical domain. 

By fixing the attention on a set S, let us express p and 
A in terms of two new variables 



p = P{u,v) X = A{u,v) 



(48) 
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We will assume that the transformation defined in (1481) 
is canonical, i. e. it has the property that there exists a 
function F{u, v) such that 



PdA = udu + dF 



(49) 



By considering u and v as fields, u = u{q,t) and v — 
v{q, t), the system (P^ . (H5|) will be transformed in a new 
system involving u = u{q,t) and v = v{q,t), having still 
a canonical structure. 

A necessary and sufficient condition that (|48l) be 
canonical is 



d{u, v) 



dP dA 
du dv 



dP dA 
dv du 



-1 



(50) 



Now, suppose that pd^(p) is singular at p = 0. Since we 
consider pj and jd, which appear in the global balance 
equation as global field variables, we deduce from 

(HT|) that we must have in any case that pd{p) is regular 
at p = 0. This means that pd'^{p) must have a pole of 
first order at p = 0. So we can write 



{a/2f 



+ g{p) (a>0). 



(51) 



where a is a constant and g{p) is a smooth function of p, 
regular at p = 0. 

By considering the canonical transformation (j48p . we 
have 

—pdtX—He = —u{q,t)dtv{q,t)—T-l'e{u,du,v,dv)+dtF{u,v) 

(52) 

where, according to ((5T|) . the transformed effective 
Hamiltonian density is given by 



n' 



—m 
2 



(a/2)2 



- 2 



( dv 



p 

{a/2)l(dP 
dv 



dP\ 
'd^J 



dkudiu 



P 

dA\ 

dv J 



-m''g{P)dkPdiP 



P 
PV . 



dkvdiv 

dP 
du 



dP 

dv 



dkudiv 



/|Ay^(a^/apy ^ /OA 

\du J P \ du J \dv 

a 

2 ■ 



(a/2)2 zap' 
P [dv , 



(55) 



So the singular dependence on P, at P = 0, disappears. 
Now, let us determine the above canonical transforma- 
tion. From (154]) it follows that logP and A/ (a/2) must 
be the real and the imaginary part of an analytic function 

logP+^^=fiz), (56) 
a/ 2 

with z = u + iv. Furthermore, due to (|50p . we have 

2 / o A \ 2 



d A 
du a/2 



1 



P-a/2 



(57) 



Then, by making use of we deduce that f{z) satisfies 
the equation 



dz 



1 

2^ 



Therefore 



-) e^^(z + ci + «:2). 



(58) 



(59) 



' where ip,ci,C2 are real constant. 

In conclusion, the canonical transformation which sat- 
isfies ((5l)) is given by 



P(u,v) 



ju + Cl)'^ + {V + 02^ 

2a 



A{u, v) 



arg{u + ci + i{v + C2)) -|- (p 



(60) 



(53) 



Our aim is to eliminate, through an appropriate choice 
of the canonical transformation, the singular dependence 
of We on P, at P = 0. If this happens, the field variables 
involved in this transformation can be promoted straight- 
forwardly to global variables. To this end, on the basis 
of the structure of the coefficients of the quadratic form 
appearing in ([53]) and of the condition ([50]) . we consider 
a transformation such that 



(In E we can select a determination of arg). So, as a 
consequence of (1511) . in this case the probability density 
is related to the squared modulus of a complex ampli- 
tude, while A is connected with the phase of this am- 
plitude. Such a situation will be called in general the 
extended quantum domain. The constant ci and C2 
can be absorbed in a further canonical transformation: 
u' = u + ci , v' = v + 02- Furthermore it is useful to use 
the scaled field variables 



adP_ _ pdA 
2 du dv 



■01 = i 

i2a)l 

so that ([50)1 becomes 

P(Vl,V^2)-V^?+^2', 



(2a)^ 

A(Vl>2) 



arctan 



If (l50l) and (l54l) are satisfied, we have 



1p2 



p(dA\fdA 
\du J \dv 



(«/2)^ 
P 



dP 
du 



dP 
dv 



= 



By considering the effective lagrangian 
Ce = -pdtX - He , 



(61) 
(62) 
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which gives the system (|45|) . (|45|) . we deduce, from 
(I521),(I531),(I511),(I5S1) and that is given also by 

Ce= - a{ipidtip2 - ip2dtipi) + 

- V{q){i;l+ijl)- 

- 2g{'4,l + VDm'^' (^i9fc^i + 4'2dk4'2) ■ 
■{ipidiipi + il>2diil)2) , 

(63) 

where we have neglected a term involving a total time 
derivative. 

Although we have derived ([63)1 in a set S, we see 
that all the coefficients of the quadratic form involving 
dk^i,di^p2 are regular functions, even at p — ipf + ip2 — ^■ 
So there is no effective restriction to consider Ce , as given 
by ((63)) . on the whole M" x [0, r]and to promote straight- 
forwardly il)i{q,t) and ip2{q,t) to global field variables. 
From (|63p we can deduce a system of two equations for 
■01 and ip2, which globalizes the system (I44|) . p5)) . 

In this respect, it is useful to consider a further canoni- 
cal transformation, by introducing a complex field ■0 and 
its conjugate 

= 01 -I- i'02 ■0* = 01 - j0'2 (64) 

In terms of the field variables and 0^* , the Lagrangian 
(|63| becomes 

Ce= + z|(0*5t7/.-Mr) - 

- ia2TO'='afc0*a,0 - 1/0V 

(65) 

From (j65p we deduce a system of two equations for tJj and 
"0*. One of these equations is the complex conjugate of 
the other. In general, there is a coupling between the two 
equations. The usefulness of the canonical fields and 
0*, stems from the fact that they allow to determine a 
single and relevant circumstance in which we obtain two 
decoupled equations. This happens when 

ff(0*0) = (66) 

In this case, we have a situation analogous (but more 
symmetrical) to that of the classical domain, in which a 
canonical field is decoupled from the other. In the classi- 
cal domain the decoupling is responsible for the existence 
of deterministic trajectories. The decoupling obtained 
when ((66)) is satisfied, is responsible for the linearity of 
the resulting equation for tjj (or 0*). In this case we 
speak of the quantum domain. When ([55]) is satisfied, 
we deduce from (|65p 

ladti^iq^t) = -^d,(m'\q)dMq,t)^ + V{q)^b{q,t) 

(67) 



For 

a = h, (68) 

(|67p is the Schrodinger equation which governs the quan- 
tum mechanics of a natural mechanical system. If 
f/(0*0) ^ 0, we obtain a non linear Schrodinger equa- 
tion. 

It is useful to stress that the need for the field variables 
like tpi{q,t) and 02(g,t) is related to the fact that our 
problem is to determine the global field variables p{q, t) 
and j{q, t), while at an intermediate step we have to intro- 
duce a local (in principle) field variable, that is X{q,t). 
We are then faced with the solution of some matching 
conditions, which among other things would require a 
knowledge (the set N) not given a priori. These condi- 
tions are automatically solved by considering new field 
variables for which there is no obstacle to be extended 
globally in straight-forward way. Without any condition 
we could not say that system of local equations (M)) and 
(l45l) is equivalent, for g{ip*ip) = 0, to the Schrodinger 
equation, which is valid globally. As a matter of fact ([33]) 
and (|i5l) is a system of generalized Madelung hydrody- 
namic equations (isl . [isj . The inequivalence between the 
Schrodinger equation and the Madelung equations has 
been already discussed by T.C. Wallstrom [2^]. However 
it is significant that, when (I66p is satisfied, the system 
of local equations (|44| and (|45| can be interpreted in 
terms of classical Brownian trajectories, as shown in the 
Nelson's stochastic mechanics approach [l^ [13, [H, [13] ■ 

Our previous analysis provides, as a particular rele- 
vant case, a scheme for a quantization procedure. In this 
respect we note that, in the framework of the canonical 
quantization rules, we have in general a problem with the 
ordering of the product of non commuting operators, as 
happens when an Hamiltonian of the type (|39p is gener- 
ally considered. Our scheme, through the equation (|67p . 
allows to give an answer to this problem. Furthermore we 
remark that the consideration of the Lagrangian ((ESI) 
or (|65p ) is also useful from another point of view. By ap- 
plying to Ce the Noether theorem Q, we can select the 
relevant conserved physical quantities, associated to the 
symmetry properties of our mechanical system. When 
(1661) is satisfied, these quantities are in correspondence 
with the expectation values of the operators of the stan- 
dard quantum theory (Tsj . 

From the theoretical point of view, our approach pre- 
dicts different dynamical behaviours, according to the 
structure of d{p). This structure cannot be given a pri- 
ori, without further considerations. In principle it must 
be argued on physical grounds. 

The quantization procedure afforded by the previous 
scheme requires that the mechanical system be natural, 
a property not necessarily contemplated by the canonical 
quantization rules. We could also develope a "quantum" 
theory of a non natural mechanical system, on the ba- 
sis of the action (l?T|) or the system ([Ml) . (|55)) . The re- 
sulting "quantum" dynamical behavior would break the 
time reversal invariance of the classical level. We note 
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that a mechanical system with relativistic kinematics is 
non natural. In this case we are also faced with an "ho- 
mogeneous" problem in the calculus of variations 0, [l^ . 

In the case of "homogeneous" problems, or problems 
in parametric form, the condition (|16p is not satisfied. In 
order to treat these cases (singular lagrangians), within 
our approach, we need a further elaboration as done in 
the calculus of variations or in the Dirac approach to 
constrained hamiltonian systems [29j . The limiting pro- 
cedure of H. J. Rothe [201 could be also useful. These 
aspects will be investigated elsewhere. 

We conclude this section with a provisional discussion 
of another extension of the classical scheme, which can 
be formulated in the framework of our approach. 

Suppose to consider the time evolution in the inter- 
val [0, r] of a dynamical system whose configurations are 
described by a discrete variable assuming a finite num- 
ber N of values. These values will be labelled by a 
discrete index running from one to N. At each instant 
t G [0, r] we consider the configuration of the system a 
discrete random variable described by a set of probabili- 
ties (a = 1,2, ..,iV). 

As a first step, in order to take into account of the fact 
that 



N 



does not depend on we introduce a quite general bal- 
ance equation in the form 

N 

+ =0 (a- l,2,..,iV) 

13=1 

(69) 

The non local and discrete equation ([Mjl takes the place 
of (fTO|) . We plan to determine variationally Pa{t) and 
la(3{t)- In this respect, here we will assume that the 
functions Pait), for each a, are such that they have at 
most a finite number of zeros in the interval [0,t]. As a 
second step we introduce a "Lagrangian" C{^), a smooth 
function of a real variable ^, and a real symmetric matrix 
Uq/3, time independent, which will be connected with the 
extension of the local substitution We will assume 

that 



and 



ci-0 = m 



Uq/3 = U/3a 7^ Va,/3 



(70) 



(71) 



Suppose now that Pa (t) and give an actual temporal 
evolution of our system and let us consider the family 
of subset r of [0, r], formed by finite unions of disjoint 
closed internals in which Pa(t) > 0, for every a. As a 
final step, in each F we extend pT|) by the "non local" 
substitution 



Pimc.l3pj{t) 



(72) 



Then we replace the action functional ([12)) by 

7a/3(0 



TV , TV 



"T a=l'^B=l ^Pi 



pima^ppjit) - 

Ac {t) ( Pc (t) + J2 (7a/3 (t) - IPa (t)))] (73) 



/3=1 ^_Pa (i)UQ/3P^ (t) 
N 



13=1 

By applying straight-forwardly the generalized action 
principle we deduce from ()73p 

Ac(t)-A^(t)^^f , ) (a,/3 = l,..7V;ter) 

\pim^f,pj{t)J 

(74) 

By introducing again a Legendre transformation 

H{v) = - m , (75) 

we can write 

dn 



-fo.f3{t) = pI{t)\Jo.f3pj{t) — {K{t) - \p{t)) 

{a,l3 = l,..N;teT) (76) 
Then we obtain the system of equations: 

A„-5IUa/34^(A„-A^) = (77) 

/3=1 Pa 

Pa + 2 ^ pi {t)\J^f3pl (t) — (A„ - A/3) = 

13=1 

(a,/?=l,..iV;tGr), (78) 

where we have made use of ([70)) . ([7T]) and ([75[) . 

The system (|77[) . ([78[) defines an Hamiltonian structure. 
It can be deduced directly from the action functional 



Ar[\a,Pa\ 



N 



a=l 



J2 PS.it)Va^pjit)n{Xa- Xf3) 
a, (3=1 



(79) 



We see that the local character of the system (1771) . ([75[) 
is determined by the non-analytic square root depen- 
dence on Pa{t) of the action functional ([79[) . Then we 
can obtain an explicit description of a global dynamical 
behavior, if we can absorb this dependence in a canon- 
ical transformation involving variables which can be ex- 
tended globally without any condition. However, in order 
to carry out this plan, we have to select appropriately the 
structure of TLir}). As a matter of fact our goal can be 
realized quite easily, by choosing as our "Hamiltonian" 



7i(r/) = h cos 



77 



(a > 0) 



(80) 
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where a and b are real constants. We have then imme- 
diately the appropriate canonical transformation, which, 
written in the standard form involving complex quanti- 
ties, is given by 



The action functional (l79l) becomes 



(81) 



r . — 1 



N 



a, 13=1 



(82) 



In terms of the variables il^ait) we obtain straight- 
forwardly a global description of the dynamical behavior. 
As a matter of fact the 4'a{t) which are solutions of the 
evolution equations are analytic functions of F, so that 
our assumption on the zeros of Pa{t) is satisfied. 

We consider two generalizations of the previous result. 
We can introduce a possible constant shift in the sta- 
tionary point of 7Q.^(t), by adding in the integrand of 
equation (I73p a term linear in ^apii) 



N 



(83) 



a, 13=1 



where the given real matrix Oap is antisymmetric {9af3 = 
—9f3a). On the other hand, through an appropriate lim- 
iting procedure, we can consider also the case in which 
some elements of the matrix Uai3 are equal to zero. Then 
we obtain the final global effective action functional 



N 
a=l 



N 

E 

a, 13=1 



(84) 



with a generic hermitian matrix given by 

Kf3 = -h\Jai3e ' » (85) 
we deduce the standard discrete Schrodinger 



From 
equation 



N 



ih^a^it) (a = 1,2, ..,7V), (86) 

p=i 

which governs the quantum dynamics of a discrete ran- 
dom variable, like the spin. 

We note finally that it possible to consider also the 
situation in which we have both continuous and discrete 
random variables. 



IV. PROBABILISTIC VARIATIONAL 
APPROACH TO CLASSICAL FIELD THEORY 

Let us come back to our physical system described clas- 
sically by a set of real fields q^{x) and the Lagrangian den- 
sity By following the theoretical scheme of SectlTTl 



we adopt at the start a probabilistic point of view and, at 
the same time, we plan to generalize the standard varia- 
tional principle based on the action ([2]). 

At first sight, the immersion of the deterministic classi- 
cal field theory in a probabilistic scheme could be realized 
immediately if a field is considered as a m.s.i. As a mat- 
ter of fact the mechanical transcription allows to obtain 
an Hamilton- Jacobi equation. In this case, we have a 
functional equation, involving a functional S(^[q{x),t]), 
for each t, of the configurations q{x) of our system of 
fields (x = {x^ , x^)) . However, we meet with diffi- 
culties if we want to place at the side of S{[q{x),t]), a 
functional probability density p([q{x),tfj . There is no 
Lebesgue measure in an infinite dimensional space and 
there are problems with a functional divergence theorem. 
So we are faced with a probabilistic scheme which is al- 
ready ill-defined at the classical level, before the tran- 
sition to the quantum theory. Furthermore, space and 
time are treated in an asymmetric way. On the other 
hand, as will be discussed in the following, it is possible 
to embed deterministic classical field theory in a simple 
and well-defined probabilistic scheme in which the above 
difficulties are avoided. 

At each space-time point x, we treat the values q = 
{q^,q^, .. , g") of the set of our fields as random variables 
described by a probability density p{q, x). We have now 
that the total probability 



p{q,x)dq 



does not depend on the space-time point x. We assume 
again that this independence can be expressed in a local 
form 

d^p{q,x)+dk3l{q.x)^Q (m- 0,1, 2, 3) (87) 

So, together with p{q,x), we have now four probability 
current densities jf_,{q,x) ={jj,{q,x), .. J^{q,x)) through 
which the " state" of the system is characterized. In terms 
of p{q, x) and jfj.{q, x) we can extend in a natural way to 
fields the generalized Hamilton's principle of scctlTTl 

It is formally convenient to take into account of the 
property that the Lagrangian density ([T]) is a function of 



q and of An independent variables — {wj^, . 

C ^ C{q,w^) 



(88) 



Now, suppose that p{q, x) and jp,{q, x) describe an actual 
evolution of our system in the space-time. We consider 
the family of subsets 17+ = U x Z, where U and Z are 
compact domains of M" and M, respectively, such that 
p{q,x) > for q,x G In each 12+, we make the 

substitution (analogous to (fTTj) ') 



Wn 



p{q, x) 



in the given Lagrangian (j88p and then we consider the 



11 



action functional 



piq,x) 



+ y{d,p{q,x)+d^jl{q,x)) 



(89) 



where, in order to take into account of the equations 
([87)) . we have introduced four Lagrange muhipher fields 
X'^{q,x), defined in Now we demand that the ac- 

tual p{q,x) and j^{q,x) are such that /n+ [jVi, P, A'^] is 
stationary under smooth independent variations Sp{q, x) 
of p{q,x) and 5i^{q,x) oi j^{q,x), with Sp{q,x) = for 
q,x e dU X dZ and 6j^{q,x) — for q e dU. Further- 
more the equations ([57]) must hold. 

By considering the variation of jf^ {q, x) we obtain from 



dkX''{q,x) 



dC 



dwj^ y p{q,x) 

(fc = l,2,..,n; i^ = 0,l,2,3) (90) 

The variation of p gives 



Now let us assume that 
act 



piq,x) dw^ \ piq,x) 



p{q, x) 







^0 



(91) 



(92) 



Then we can solve (1901) with respect to jfi/ p. Through 
this solution the equations (j87p and ((9T|) can be expressed 
in an appropriate and useful form, by making use again 
of the Legendre transformation. 

Let us introduce a generalized " Hamiltonian" function 

n 



(9J 

where 'w^{q, tt^) is the solution of the system of equations 
„ dC 



dwl 



(93) 
tions 

(94) 



Then we have 
P dK 



{q,dXn (aA'' = (9lA^..,9„A'') (95) 



So we obtain the system of equations 

d^X''{q,x) + n{q,dX''{q,x)) = (96) 

/ d'l~L \ 
df,p{q,x) + dk(p-^{q,dX''{q,x))j^O 

{q,xen+) (97) 



Due to the absence of an explicit restriction connected 
with the set of points in M" x M where p{q, x) = 0, we 
can reconsider ([96|) and ([97|) as equations in R" x M, free 
from constraints. They are the version, in the framework 
of the previous probabilistic scheme for classical field the- 
ory, of the equations (PT|) and for a m.s.f. Through 
([M)) we have obtained a generalization to field theory of 
the Hamilton-Jacobi equation. However this generaliza- 
tion is again a partial differential equation, and not a 
functional equation, as happens in the standard canon- 
ical field theory. Furthermore the space coordinates are 
treated on the same footing of the time. These four pa- 
rameters together take the place that the time alone had 
in previous sections. The system ([M|) . ([571) satisfies man- 
ifestly relativistic invariance, if C, as we will always as- 
sume, is a Lorentz scalar. The generalized "Hamiltonian" 
7i, given by (|93p . is not an energy density, but a Lorentz 
scalar. Another feature of our approach to field theory, 
is the occurrence of a system of four equations for p, in 
place of the only one equation for a m.s.f. As a 
consequence (^5)) is now not decoupled from ([W)) . The 
functions X^{q,x), besides to satisfy (IM)) . have to guar- 
antee the consistency of the system (1^71) . 

The Euler-Lagrange equations of the standard deter- 
ministic classical field theory, can be considered as a par- 
ticular consequence of the system . ([77)1 . 

Let q{x) — {q^{x), .. , g"(a;))) be a classical field. It can 
be easily verified that ((97)) admits as a particular solution 
the probability density 

p{q,x)=Siq~qix)) (98) 

if, and only if, the field q{x) is such that 



d.qHx)^^[q{x),dX^{q{x),x) 



(99) 



On the other hand, by deriving (I96p with respect to q^, 
we see that q(x) must satisfy also the equation 

d.d.X{qix),x) = -^(^q(x),dX^iq{x),x)^ (100) 

By taking into account of the Legendre transformation 
dMl), we deduce from ((99|) 



d,X^iqix),x) 



dC 

dw], 



[q{x),d,q{x))=^'i{x) (101) 



Then the equation (jlOO[) gives 

that is the Euler-Lagrange equations ^ for the field q{x). 

A function q{x) which is a solution of ([M]) . with 
A''(g,a;) solution of (196]) . is said to be embedded in a 
geodesic field d, @, [25]. As we have seen, such a func- 
tion is an extremal, i. e. a solution of ([3]). It can 
be shown that any extremal q{x) can be embedded in 
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a geodesic field, in some neighborhood of the surface 
q = q{x) UMMMM- 

We can take directly as unknown the fields 
diX^{q{x),x) = 7r['(a;). Then we can consider directly 
the system 



dn 



q{x),Ti''{x) 



d,n';{x) = -^(^q{x),7r^x) 



(102) 



(103) 



which is the rclativistic covariant version, in classical field 
theory, of the canonical Hamilton's equations. 

As a matter of fact, the equations . . 
(I102p . (|103p . which we have deduced through our prob- 
abilistic approach, have already been obtained on the 
basis of different considerations, in the deterministic 
De Bonder- Weyl theory developed for multi-dimensional 
variational problems Q, 0, i, H H, [ll HI ■ In the frame- 
work of the deterministic classical field theory, they have 
been also reobtained elsewhere [l^, [l^] • 

We note that, in the De Bonder- Weyl theory for clas- 
sical fields, the relevant aspect of manifest relativistic 
covariance compels to have four "momentum density" 
TT^{x) conjugate to each classical field q'^{x). This makes 
problematic the standard Poisson brackets formalism and 
then the transition to a quantum theory. However this 
redundance of conjugate variables 7rf (x)can be reduced 
through proper constraints. As a matter of fact these 
constraints are contained in the equations (|102p . 

We can consider a family of space-like surfaces which 
covers simply the manifold M. Such a family can be 
related to a description of the dynamical evolution of our 
system. However, in the following, for simplicity, we will 
limit ourselves to surfaces having xq constant. 

Now the three equations in p02p 



dniq\x) 



dn 



g(x),^'^(x),7r"(x)) (m = l,2,3) 

(104) 

can be interpreted as constraints for the fields 
7r"(a;) (n = l,2,3). 

By solving them, we will obtain 



i:T{x)^Fr[q{x)y{x\d^q{x) 



(105) 



Then we can introduce the function 
'Hc\q{x),-n^{x),d^q(x)\ = 



= n[q{x), 7t"{x), F"'{q{x),7r'{x), d„q{x)) 

-{d^q\x))Fr(^q{x),7r^{x),dnq{x)'j (106) 

We have 

Fr{q{x),7:^{x),dnq{x)) 



did^mq^x)) 



(107) 



The system (|T02l) . (|T03l) is then reduced to 



doq\x) 



do7r°{ 



fdUc 



\ dqi 



d„ 



d{d^q'{x)) 



(108) 



(109) 



that is the system of equations of the standard Hamilto- 
nian formalism Q , related to the description of a classical 
field as a m.s.i. Really, it follows from ([M]) that Tic is the 
standard Hamiltonian density. 

Then we see that, if we carry out the elimination of 
the conjugate fields 7r"(a;), the resulting theory (which is 
the standard one), meets with the difficulties mentioned 
at the beginning of this section. On the other hand, if we 
keep the complete manifest covariant structure of the Be 
Bonder- Weyl theory, then this theory, as we have shown, 
can be embedded in a simple and well defined probabilis- 
tic scheme, related to a relativistic invariant generalized 
partial differential equation of Hamilton-Jacobi type. In 
such a scheme it is possible to contemplate quantum ex- 
tensions, as in the previous section, which arc disengaged 
from the usual s.c.q.r. This will be done in the next sec- 
tion. 

Bue to ([59)1 and ([55]) , the basic equations of the previ- 
ous scheme, that is (jHH) and ([57)l . can be derived directly 
from the standard variational principle applied to the ac- 
tion integral 



Ao[p,A^] - / dxdq 
Jn 



-y{q, x)d,p{q, x)+n{q, dX^piq, x) 



(110) 

where ft is now a generic domain in of M" x M. The 
variations of must vanish on the boundary of f2. From 
the point of view of standard classical field theory, we see 
then that only X'^{q,x) has a conjugate variable, that is 
p{q,x), while for X™{q,x) (m = l,2,3)such conjugate 
variables are absent. Therefore three of the equations 
(l97l). that is 



d,Mq,x)+dk(^p^{q,dX''{q,x))^ =0 (m = l,2,3) 

(111) 

have to be considered as constraints of the dynamical 
evolution, in analogy with what happens for the system 
(|102p . We note that the constraints which appear in the 
previous approach are intrinsically related to the prop- 
erty that the generalized hamiltonian is a scalar, and not 
an energy density. As a matter of fact, due to ([5^ . we 
are not treating singular Lagrangians. It seems that these 
constraints have a role different from that given by the 
Birac approach. 

We conclude this section with the following remarks. 
First of all, the previous probabilistic approach to classi- 
cal field theory requires, in general, some caution. This 
is related to the condition (|92p . which makes possible a 
transition from the Lagrangian formalism to a covariant 
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Hamiltonian formalism. In the case of a classical Dirac 
field, it is well known the Lagrangian density is linear in 
the derivatives of the field components so that the condi- 
tion ([^^ cannot be satisfied. Another problem, of tech- 
nical character, appears if we have constraints between 
the components of classical fields, as happens in the case 
of vector fields. The role and meaning of these cases in 
our probabilistic approach will be investigated elsewhere. 
In the following we will limit our considerations to scalar 
fields having the property to satisfy the condition ([9^ . 

Since we assumed that our classical Lagrangian density 
(H]) does not depend explicitly on a;, we can consider the 
invariance of our probabilistic theory under space-time 
translations, that is under the transformation 



y = {x, q) 



y 



ix + s,q) £ = (£",e\£^£-^) (112 



Then, from ()112p and a straight-forward application of 
the Noether's theorem to the action (jllOp . we deduce 
the local conservation law 



d^{T!;{q,x)p{q,x)) + d,Ql{q,x) ^ , 



(113) 



with 



T^iq,x) = dk \^ {q,x)^{q, d\^) + g:(n- d^X" |^ 

(114) 

Ql, is a function of q,p{q,x),dX^{q,x). 

So, our approach contemplates a random energy- 
momentum tensor Tj^{q^ x) given by (jll4p . If we consider 
the mean of Tlf{q, x)at the point x 

Tt{q,x)^ I Tl^{q,x)p{q,x)dq (115) 



we deduce from (|113p 



(116) 



In the case of the solution (jOS]) , we obtain from (|114p and 

(ESI) 



Tl[x) = T,^(x) 
dC 



dwl, 



q{x),d^q{x))d^q\x) ~ gZC{q{x),d^,q{x)) , 

(117) 



i. e. the standard result of the deterministic classical 
field theory. 



V. THE QUANTUM EXTENSION FOR 
SCALAR FIELDS 



We develop a quantum theory for scalar fields, starting 
from the probabilistic scheme of the previous section and 
following the approach of the sect lllll We will limit our 



considerations to a set of n real fields described by a 
"natural" Lagrangian function 



V{q) 



(118) 



where rjij (q) is positive definite Vq {rjij (q) = rjji (q) , the 
inverse matrix is denoted rj^^ (q)) .The Hamiltonian func- 
tion in this case becomes 



nq,nn = ^V'- 



V{q) 



(119) 



As a first step of an extension of the scheme of the previ- 
ous section, we take into account, in the balance equation 
([87|) . of an additional contribution to d^p{q,x). This is 
given by current densities i^{q, p{q, x) , dip{q, x) , ...) hav- 
ing a local dependence on p and on its derivatives with 
respect to q*. So the equation ((87)) is extended to 



d^p{q,x) +dkj^{q,x) + dkil{q,p{q,x),djp{q,x),..) = 

(120) 

Through we introduce dynamical mechanisms of dif- 
fusive type. These mechanisms are internal so that, due 
to the homogeneity of space-time, they are assumed to 
have no explicit dependence on x. However, due to the 
structure of ijj, we are faced also with the problem to as- 
sure that, in the theory which we are developing, the 
isotropy of space-time is not destroyed. To this end, 
we note that, as in the case of a m.s.f., the transition 
from (|87p to (|120p makes relevant the choice of a La- 
grangian density which, at the classical level, is equiva- 
lent to £[q, jij_{q, x) / p{q, x)) . As a matter of fact, when 
(|87| is satisfied, the Lagrangian density 

C = c(q,i^]+D{p)a-'^-^^dkp{q,x), (121) 



P 



p{q,x) 



where D(p) is a given function of p and a'^ is a given 
constant four vector, is equivalent to C. 

Then, as a basis for a quantum theory of our system, 
we consider the action functional 



dxdq 



C{q,^-^)+D{p)a-^-^dkP 
P P 



+ y{d,p + dkjt + dkil{q, p, dip, ..)) I 



(122) 



The equations p2ip and (|122p are the generalization 
to our fields of the equations ((26|) and ([31]) for mechanical 
systems. 

In the case of non relativistic mechanical systems we 
were faced with the breaking of time-reversal invariance, 
which is related to the choice of an arrow of the time 
axis. In the case of fields, the introduction of a con- 
stant four vector a'', that is a privileged direction in the 
Minkowski space, breaks the isotropy of the space-time. 
At the classical level we have no problem, since the ad- 
ditional term in the equation (|12ip is ineffective. This is 
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not the case, in general, if we start from the action func- 
tional (|122p . However the isotropy of the space-time can 
be still preserved by choosing appropriately the current 
i^. As a matter of fact, the request for the space-time 
isotropy allows to determine completely i^, starting from 

the knowledge of C In fact, on the basis of this request, 
we deduce from pl8)) . (fm]) and p22)) that 

it{q,p,dip,..) = a,fi''Hq)pDip)d,p, (123) 



which is a generalization of (|4ip . 

Then we deduce from piSp and (jl22p the Lorentz in- 
variant system of equations 

OA'' + ]^ifd,Kdjy -f v{q) + 

+a.a'^[]^i^^={pD^{p))'d,pd,p-d.{i^^=pD^{p)d,p)] = 

(124) 

d^p + dk{pv'''d,\^)^0 ((/i = 0,l,2,3), ((z,a;er!+)), 

(125) 

which appears as a straight-forward generalization to 
fields of the system (j44 |) .(|45 p . obtained in the quantum 
theory of a m.s.f. 

The system (|124p . (|125p can be derived directly from 
the action integral 



dxdq 



~y{q,x)d,p{q,x)+n{q,d\^')p 



-pD\p)f^^'dkpdip 



where Ji is given by (|119p (the variations of A'' must 
vanish on the boundary of 0+). 

The last term in the equation ()126p has the same struc- 
ture which we have found in the case of a m.s.f. So we 
can speak for fields of an extended classical domain or of 
a quantum domain, according to the behavior of pD^{p). 
We will fix our attention on the quantum domain. This 
will be characterized, as in the case of mechanical systems 
{eq.l^,g{p) = O) , by requiring 



avo!" pD'^{p) 



(//2)^ 



(/>0): 



(127) 



The system (|125l) . (|128l) has a Lorentz invariant local hy- 
drodynamic form, which generalizes the Madelung hy- 
drodynamic equations related to non relativistic quan- 
tum mechanics. The system is also invariant under the 
transformations involving space-time inversion 



9' 

p{q,xn 
X^{q,xn 



q'^q' 

p'{q,x"')=p{q,xn 
\'^{q,x"') ^ -X^'iq.x") 



(129) 



We note that, while in the case of mechanical systems 
the constant a of ea. ([5T|l . analogous to /, has the dimen- 
sion of an action (as a matter of fact a = h, according 
to (pB]) ). in the case of our fields the constant / has the 
dimension of a (spatial) density of an action, as a conse- 
quence of our previous equations (|126p . (|127p . This means 
that, by considering h, f can be related to a fundamental 
quantity having the dimension of a spatial length. 

We will see some consequences of the system 
(|125p . (|128p . by assuming that, as in some usual models 
of field theory (Klein-Gordon, Aq"*, ...), 



Viq) > , V{q) -l-oo for\q\ +c 



(130) 



Furthermore we will consider the random energy- 
momentum tensor Tq'^^q^x) in the quantum case. This 
can be obtained by adding to the classical expression 
(|114p the "quantum potential" which appears in the 
equation (|128p . Taking into account of pi9p . we have 



(126) TQl{q,x) = 7^^Hq)d,X''d,X.+g:(^- 



p dM^d.p^/^) 



,1/2 



+ v{q)--v'=d,x^d,x^ 



(131) 



for (g, x) e f2+. 

In particular, the random energy density e'^{q,x) 

and the random momentum densities P^{q,x) (m = 
1, 2, 3), are given by 



eQiq,x)^TQ°{q,x) 



i7y'^9,A"9,A" + ir7*^9,A"a,A™ + 
V{q)- Y 1/2 (132) 



where / is constant. We have assumed in (|127p that 
a^a" > 0, i. e. the currents are time-like, as can be 
expected on physical grounds. 

In the quantum domain (|124p becomes 



d.x'' + -Ti'^ d.x^d.x" + v{q) + 



pfi ,., d,pd,p 1 



4 V2 



d^iv^'d.p] 



= d^X" + -Tj'^A.djX" + V{q) - 



pdM^djp^/^) 



0, 



(128) 



P,^(g,x) = rQ°„(g,x) = 77^^-9feA09,A™ (m = 1,2,3) 

(133) 

The simplest solutions of the system (|125p . (|128p can 
be obtained by considering the case 

a,A^(g,x)=0 (i= l,2,..,n; /.t = 0,l,2,3) (134) 



In such a case, (|125p is solved automatically by p{q, x) 
independent of x, that is p{q,x) — p{q). Furthermore 
(I128P (with p) requires necessarily that also dtyX"^ is in- 
dependent of X, that is, due to (|134p . 



d^X^iqjX) = -w , 



(135) 



15 



where w is a constant. Then (|128p becomes 

+ V{q)^w (136) 



21/2 



i 

for q such that p{q) > 0. 

This equation can be translated in a global form, by 
setting 

(137) 

so that we are led to replace (|136p with 



- y5.('7'''5j-V'('?)) ) + ViqWq) = w^{q) , (138) 



that is a stationary Schrodinger equation with the con- 
stant /, instead of h 

Let us denote by Wr the eigenvalues of the " Hamilto- 
nian" operator 



Z! 

2 



H = -^d,W^dr) + Viq). 



We order them according to 



< Wo < wi < 



We conclude that our system (|125[) , (|128[) admits as so- 
lutions invariant states, which are associated to the eigen- 
values Wr of H. In an invariant state, corresponding to 
a particular eigenvalue Ws , the energy momentum tensor 
is well defined, constant and given by 



(139) 



Due to (|139p . we see that it is a consequence of our 
equations the existence of several non-degenerate vacuum 
states, each having a constant finite energy density. We 
will call the state corresponding to wq the fundamental 
vacuum, while the other states associated to Wr [r > 0) 
may be called virtual vacua. The energy gap between the 
fundamental vacuum and a virtual vacuum is infinite. 
We note that, in our approach, there is no substantial 
distinction between the case of interacting fields and the 
non interacting one. The fundamental vacuum exists also 
for interacting fields. In the traditional approach, in the 
case of interacting fields, the existence of the vacuum 
state must be, in some way, postulated. In the case of 
free fields, the traditional approach gives an infinite zero 
point energy, but this is infinite also in any finite spatial 
volume, while our wq is finite. 

The previous aspects of our approach may be relevant 
in connection with the dark energy problem [28| . We 
note also that, as a consequence of (|130l) . the eigenfunc- 
tions ipr{q) associated to the invariant states are rapidly 
decreasing for \q\ — > +oo. Then we have that the fluctu- 
ations 

(Qi - Qi)iqj - Qj) = j dq{q^ - qi){qj - qj)ifl{q) (140) 



{qi = J dqqiipQ{q), the ipr{q) are supposed normalized) 
are finite. In the traditional approach, if we consider, 
for example, the real Klein-Gordon field, the analogous 
quantity {0\q^{x) |0) (q the operator valued Heisenberg 
field) is infinite. 

To simplify our further discussion, in the following we 
will limit ourselves to the case of only one scalar field, 
assuming also that riij{q) = rj{q) = r/ (positive constant). 

Now, as next steps, we will examine space-independent 
solutions and time- independent solutions. The first type 
corresponds to solutions such that 



d_ 

dq 



A"((?,a;) = (m = l,2,3) (141) 



but / dq 7^ 0, in general. 

With p{q,x) — p{q,x°), as a consequence of (|14ip . the 
system P^ . p^ becomes 



a|rAo(.,x) + J^A'"(,,x) + i-(|Ao(,,x))V 



(143) 



From (fT42ll we deduce X°iq,x) = A(g,a;°), while (fT4T|) 
and ^43]) require that 5|^A™(g,a;) = f{x°). Writing 

we can absorb A(a;'^) within X{q,x'^), by considering 
X{q, x°) = X{q, x^) + A{x^). Then we arrive at the equa- 
tions 



2-1] \ dq 

dUSI) 

According to sect Hill this system is the local form of 
the global linear Schrodinger equation 

a;0) = -^^^{q, + V{qMq, x^) , 

where locally 



^{q,x'^)^p'^\q,x")exp{yXiq,x") 



(144) 



The general solution of (|144p is a linear superposition 
of the vacuum states. Due to (|133p and (I14ip . we have 
for this solution, P^{q,x) =0 {m = 1,2,3), while the 
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random energy density is, in general, time-dependent. 
On the other hand, its expectation value is a constant. 
In fact, due to (fT43f|l . (fT33| gives 

(145) 

so that 

dqe'^{q,x'')p{q,x°)^w = J dq^ {q, x^f [q, x'') 

(146) 

Now we discuss the very interesting case of time- 
independent solutions, by considering the situation 



dq 



{q,x) = 0, 



(147) 



while, in general, -^\"^{q,x) ^ 0. Also in this case 

P^{q,x) — 0, but now there appear important features 
of the random energy density. 

As a consequence of (jl47p we have p{q, x) = 
p{q,x),X"'{q,x) = A™(g,f), (f = {x^,x^,x^),m = 
1,2,3) while, taking into account of (|128p . we can write 
g|oA°(g,x) — g{x). By introducing h"^{x), such that 
g^/i™(x) = g{x), we can absorb h"^{x) within \"^{q, x), 

by considering X'"'(q, x)+h"^{x), which we recall X"^{q, x). 
Our system (jl25l) . (|128p becomes 

-Q-^P{q, f) + -— f /5(q, x) — Xrniq, x)j ^0 (m = 1, 2, 3) 

(148) 



This cannot be translated in a Schrodinger type equation, 
due to the minus sign in the third term of (|149[p . However 
it is useful to introduce a transformation involving two 
real, "conjugate" , functions, as is done in diffusion theory 
[2]| . We consider two real functions, </>((?, r) and ^(q, r), 
such that 



p{q,r) = 4>{q,r)(t){q,r) , 



(151) 



Mq^r) ^Iwor log iwhenp>0) (152) 

3 2 (t>{q,r) 

In terms of (t){q, r) and (i>{q, r), the system pigft . ljnWI) 
becomes 

+ ^logf4)0(,,.) (153) 
r (j){q,r)J 



d ~, 



P d 



fg^<f>{q,r) ={ -L-—+V{q)-Wo + ^\o^ 



2rj dq 



<i>{q-, r) 



4){q, r) 
(154) 

We are interested to bounded solutions of this system. 

We note that, if we neglect the term 
(//r) log((/)(g, r) / 0(g, r)) , the resulting equations 

- /|:0o(g, r) = ( - ^ ^ + V{q) «;o) Mq, r) (155) 



d 



dx 

+V{q) 



p^/2(g,x) = (149) 



2r;pi/2 dq^ 



In order to see some basic aspect of this system, we 
fix the attention on its possible spherical symmetric so- 
lutions. We start from the ansatz 



A"(q,f) = -A"(q,r): 



(r = |f|,m= 1,2,3) 



(150) 

We infer then from l|148p p{q, x) = p{q, r). The system 
P15)) . pig)) becomes 



f-^Mq, ^) = ( - + ^('?) - ^0 ) Mq, (i56) 



can be easily solved. The ea. (|155l) . which is a diffusion 
equation, admits non negative solutions, bounded for ev- 
ery r, having the structure 

hiq^r) = co^po{q) + ciV'i(g)e~' + ... 

+ c„Vn(g)e""""j^™"''' + ...(157) 



On the other hand, (|156p admits only one solution, which 
is bounded for large r, that is 



-|:A(,,.)-'a(,,0-^(|a(,,0) +viqy 



2rip^^^(q, r) 



■p^/\q,r)^Q 



(HSI) 



0o(g, ?■) = coV'o(g) (158) 

Now we take co — co 0), so that, for large r, 



e f ' 



L log = /fli^^e-'^^+O 

r Mq^r) cor^o(g) V r 

(159) 

(assuming ci 7^ 0). Due to ()159p . we can apply to the 
system (|153p . (|154p . for large r, the method of successive 
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approximations, starting from (jjo^q, r) and (t>o{q, r). As a 
first improvement we obtain 

nq,r) ^ 0o((?,r)+ / log^- -dr = 

Jr f 0o(9, 

/ 1 



coMq) + ciV'i('?) 1 + 



P 1 



,r) ~ CQ^po{q) - Ci-0i(g) 



tui — Wo 1" 
6-^^^ + ... (160) 



/ 1 



wi — wq r 



P 1 



('"1 "'"O) „ 

e 7 + 



(161) 



The last dots represent terms of order exp{—2{wi — 
wo)r/ f) or exp{—{w2 — wq)/ f), or smaller than these. 
From (fT5T|) . (fTH0l) . (fT6T]) we obtain that, for large r, 



p{q,r) 2:^ipl{q)+Ci'iJjo{q)'iJji{q)e 7 + 



(162) 



where we have taken cq = +1, for a proper normalization. 
Since ipoiq) is the state of the fundamental vacuum, (I162|) 
describes a situation in which the physical system of our 
field is static and confined within a region of the space, 
centered at the origin, of radius ~ f /{wi ~ wq). In such 
a situation we have that, besides the constant vacuum 
energy density, there is a random energy density confined 
within a sphere of radius f /{wi — wq). 

We obtain analogous results if ci = and C2, for ex- 
ample, different from zero. We note that, due to the 
invariance of our equations as regards spatial transla- 
tions, the previous considerations on static solutions can 
be applied to a generic sphere having center at an ar- 
bitrary fixed point of the space. Then one can consider 
also several spheres, with centers having distances larger 
than 2//(u'i —wq). The above results seem to have some 
appealing aspects. Their interpretation and a further 
analysis of our equations for specific models will be given 
elsewhere. We note that in the case of mechanical sys- 
tems, where we have only the parameter t, the balance 
between two competitive diffusion processes makes inef- 
fective the H-theorem. On the other hand, as we have 
seen, in the case of fields this balance does not prevent 
that an H-theorem be operative in the space-like direc- 
tions. 

We add a brief comment on the dynamical behavior of 
our system (jl25l) . (|128l) . We can also introduce for fields 
a quantity like the wave function of the mechanical sys- 
tems. However, in this case, the notion of wave function 
is associated to some covering of the space-time through 
a family of space-like surfaces. Each such a covering de- 
termines a particular wave function. To fix the ideas, 
let us consider the standard family of space-like surfaces, 
each having a constant "time" xq. To this family we can 



associate the wave function il>(q,x'^\x), parametrized by 
X, given locally by 

^{q.x'^lx) = p^/^{q,x°\x)exp{jX'^{q,x'>\x)) 

Due to sect lllll the system (I125p . (ll28p can be trans- 
lated in the system 

tfJL^{q,x'\x) = -ild,(^r^^^q)dMq.^°W?) + 

+ (^A'"(g,a;0|f) + y^\q)d,X^iq,x'^\x)d,X"^iq,x^\x) 



+ V{q)j^P{q,x°\x) = 



(163) 



-{i:ij)+dk{-4>iJv'''d,X.m)^0 (m = 1,2,3) 



dx' 

m) 

We have then for ipiQj x'^\x) a Schrodinger equation, in 
which, besides the known function V{q), there is also a 
"self potential, constrained by the equations (I164p . As 
a result, even if we introduce the above wave function, 
the linearity found for mechanical system is lost in the 
case of the fields. 

However it is possible that linearity requires a mathe- 
matical object more general than ^Jj(q,x'^\x). 

We conclude here with some few remarks. The pre- 
vious approach for scalar fields is based on a particular 
well-defined probabilistic scheme, which makes use of or- 
dinary functions. On the other hand the standard ap- 
proach of the quantum theory of fields, is based on a 
more general probabilistic scheme which is of functional 
type, but ill-defined. However the probabilistic distribu- 
tions of our approach could be considered or imposed, 
as the appropriate marginal distributions of functional 
distributions which are not well-defined. 

Our approach requires some further developments in 
order to treat the case of singular Lagrangians. An anal- 
ysis of the Dirac fields is particularly important in order 
to see how the previous probabilistic scheme for scalar 
fields must be modified. A further important problem is 
to investigate the application of our approach to the case 
of the gravitational field. 
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